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ABSTRACT 


The  combustion  of  a  porous  graphite  fiber  plate  is  analyzed 
A  transient  one-dimensional,  mathematical  model  is  utilized  to 
conduct  numerous  computer  tests.  The  model  simulates  the  ther¬ 
mal  response  of  a  porous  graphite  fiber  plate  and  it  is  based 
on  appropriate  energy  equations.  Two  of  the  energy  equations 
are  based  on  energy  balances  conducted  on  both  the  porous  struc 
ture  and  on  the  air  within  the  graphite  fibers.  The  third 
energy  equation  is  based  on  a  mass  balance  performed  on  a 
differential  volume  of  porous  medium. 

A  test  plan  is  formulated  and  computer  test  runs  are  con¬ 
ducted  to  investigate  the  effects  of  plate  thickness,  filament 
diameter,  permeability,  porosity  and  exterior  wind  velocity  on 
the  thermal  response  of  the  porous  graphite  fiber  structure. 
The  test  results  are  analyzed  and  conclusions  are  presented. 
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I.  INTRODUCTION  TO  THE  PROBLEM 


The  advent  of  advanced  fiber-reinforced  composite  materials 
has  spurred  substantial  research  efforts  across  a  broad  spec¬ 
trum.  Whole  organizations  have  been  formed  to  analyze,  design, 
and  fabricate  components  made  of  composite  materials. 

Fibrous  composites  consist  of  long  fibers  contained  within 
a  matrix  material.  Long  fibers  are  inherently  much  stiffer 
and  stronger  than  the  same  material  in  bulk  form.  Fiber  pro¬ 
perties  differ  from  those  of  the  bulk  form  due  to  the  more 
perfect  structure  of  the  fiber.  Additionally,  fewer  internal 
defects  are  present  in  fibers  than  in  bulk  material.  Strengths 
and  stiffnesses  of  a  few  selected  fiber  material  are  presented 
in  Table  I  below: 


TABLE  I:  FIBER  PROPERTIES 


FIBER 

DENSITY. o 
[lbm/ in5 ] 

TENSILE 
STRENGTH , S 
[103lb/inz] 

S/o 

[10s  in] 

TENSILE 
STIFFNESS.  1 
[104lb/ inz ] 

E  H/o 

[  10  7  in] 

Aluminum 

.097 

90 

9 

10.6 

11 

Titanium 

.170 

280 

16 

16.7 

10 

Steel 

.  282 

600 

21 

30.0 

11 

Boron 

.093 

500 

54 

60.0 

65 

Graphite 

.051 

250 

49 

37.0 

72 

Boron  and  graphite  fibers  are  frequently  selected  for  com¬ 
posites  that  will  be  utilized  as  aircraft  structural  material. 
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The  matrix  material  in  such  applications  is  frequently  a  plas¬ 
tic  such  as  epoxy  or  polyimide.  These  advanced  composites  are 
characterized  by  properties  that  are  superior  to  those  of  ma¬ 
terials  that  have  been  previously  utilized  for  aircraft  struc¬ 
tures  (i.e.  aluminum  and  titanium).  Two  maior  advantages  are 
improved  strength  and  stiffness,  when  compared  with  other 
materials  on  a  unit  weight  basis. 

Advanced  composite  materials  have  been  designed  that  are 
three  times  as  strong  as  aluminum,  yet  weigh  only  sixty  percent 
as  much  [Ref.  1] .  It  is  this  favorable  strength  to  weight 
ratio  that  makes  composite  materials  particularly  attractive 
to  the  aerospace  industry. 

Currently,  almost  every  aerospace  company  is  developing 
products  made  with  f iber-reinforced  composite  materials. 
Replacement  components  are  frequently  made  with  composite 
materials.  Examples  of  this  application  are  the  boron/ epoxy 
fuselage  section  and  horizontal  tail  on  the  General  Dynamics 
F-lll.  Several  companies  are  also  introducing  production 
components  made  from  composite  materials  on  their  new  aircraft. 
An  example  of  this  application  is  the  F/A  18  Hornet  Strike 
FigFiter  currently  being  built  by  McDonnell  Douglas.  The  F/A  18 
airframe  is  a  balance  of  conventional  materials  and  graphite 
epoxy  composite.  The  wing  skins,  trailing  edge  flaps, 
stabilators,  vertical  tails  and  pudders,  speed  brake,  and  many 
access  doors  are  made  from  composites.  This  balanced  design 
has  resulted  in  appreciable  weight  savings  and  increased 
aircraft  performance. 
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Weight  savings  are  of  critical  importance  in  military  air¬ 
craft  applications.  Weight  reductions  result  in  the  following 
desirable  features: 

(1)  Increased  payload 

(2)  Improved  operating  efficiency 

(3)  Extended  range 

It  is  these  attractive  features  that  have  created  a  great 
deal  of  interest  within  the  U.  S.  Navy.  Lighter  weight  air¬ 
craft  with  increased  payload  are  ideal  for  carrier  operations. 
Additionally,  composites  materials  have  improved  the  feasibility 
of  VSTOL  aircraft,  which  could  lead  to  smaller,  less  expensive 
aircraft  carriers. 

The  introduction  of  aircraft  structural  components  made 
from  composite  materials  and  the  potential  of  aircraft  made 
entirely  from  composites  has  spurred  a  great  deal  of  Navy 
sponsored  research  activity.  Research  is  being  conducted  to 
gain  insight  into  the  following  features  of  composite  materials: 

(1)  Thermal  properties 

(.2)  Strength  and  stiffness 

(.3)  Fatigue  life 

(4)  Combustion  characteristics 

The  combustion  characteristics  of  composite  materials  and 
their  thermal  response  in  a  fire  environment  are  of  particular 
concern  to  the  Navy.  The  potential  for  aircraft  carrier  flight 
deck  pool  fires  is  quite  high,  due  to  the  presence  of  various 
volatile  fuels  and  ordnance  on  the  flight  deck  during  flight 
operations . 
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The  response  of  metallic  materials  to  aircraft  carrier 
flight  deck  pool  fires  and  the  required  fire-fighting  tech¬ 
niques  are  well-known,  but  little  data  is  available  concerning 
the  combustion  characteristics  of  composite  materials.  In 
order  to  gain  an  insight  into  these  characteristics  an  initial 
series  of  tests  was  conducted  at  the  N'aval  Weapons  Center, 

China  Lake  [Ref.  2].  The  objective  of  these  tests  was  to  de¬ 
termine  the  response  of  graphite  -  epoxy  composite  materials  to 
carrier  flight  deck  pool  fires  and  compare  them  with  equivalent 
strength  aluminum. 

For  these  tests,  a  composite  test  specimen  was  mounted  in 
such  a  manner  as  to  simulate  the  heating  condition  of  a  complete 
wing  subjected  to  a  large  pool  fire.  The  various  test  specimens 
were  then  subjected  to  a  steady  state  thermal  environment. 
Specimen  temperature  data  were  recorded  by  thermocouple  instru¬ 
mentation  until  all  internal  reactions  had  ceased. 

The  combustion  process  was  observed  to  proceed  in  two 
distinct  phases: 

Cl)  Epoxy  combustion 

C2)  Graphite  combustion 

At  the  completion  of  the  first  stage,  the  epoxy  was  totally 
consumed,  exposing  the  porous  graphite  fiber  structure  (epoxy 
ignites  and  burns  at  a  much  lower  temperature  than  graphite)  . 

In  some  cases,  the  remaining  graphite  fibers  exhibited 
characteristics  of  sel f - sus tained  combustion  and  a  tendency 
to  smolder,  while  in  others  the  fibers  would  cool,  indicating 
an  inherent  resistance  to  fire  burn- through . 


These  initial  tests  indicated  that  the  thermal  response  o 
a  composite  plate  is  not  constant,  but  is  a  function  of  the 
various  parameters  that  characterize  the  composite  material 
(i.e.  filament  diameter,  ply  thickness,  plate  thickness),  the 
exterior  wind  velocity,  and  the  thermal  environment.  The 
number  of  parameters  involved  in  the  problem  indicated  that 
an  exhaustive  experimental  program  designed  to  further  inves¬ 
tigate  the  thermal  response  of  a  composite  plate  would  not  be 
economically  feasible.  Therefore,  a  computer  model  appeared 
to  be  a  desirable  means  for  conducting  numerous  tests. 


II.  DESCRIPTION  OF  MATHEMATICAL  MODEL 


In  order  to  conduct  a  thorough  computer  analysis  of  the 
combustion  characteristics  of  a  graphite-epoxy  composite  plate, 
an  appropriate  mathematical  model  was  required.  A  suitable 
one-dimensional  model  was  previously  formulated  by  Vatikiotis 

/ 

[Ref.  3].  The  model  was  to  simulate  the  thermal  response,  in 
one-dimension,  of  the  graphite  fiber  mat  after  epoxy  consumption. 
The  following  discussion  summarizes  the  assumptions,  governing 
equations,  boundary  conditions , and  method  of  numerical  solu¬ 
tion  utilized  in  the  model. 


A.  AVERAGE  PORE  VELOCITY 

In  formulating  the  model  it  was  assumed  that  an  exterior 
wind  velocity  would  always  be  present,  thereby  providing  air 
flow  over  one  surface  of  the  composite  plate.  This  assumption 
is  quite  realistic  when  one  considers  ther  relative  wind  velo¬ 


cities  consistently  present  on  aircraft  carrier  flight  decks. 
The  air  flow  produces  a  pressure  differential  across  the  com¬ 
posite  plate.  This  pressure  differential  induces  a  convective 
air  flow  through  the  porous  graphite  structure  governed  by 


Darcy's  Law, 

-m  dP 


(II. 1) 


where  u^  is  pore  velocity,  y  is  the  dynamic  viscosity,  '  dx 
is  the  pressure  gradient,  and  m  is  the  specific  permeability 


of  the  medium. 
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As  a  result  of  extensive  experimental  work,  it  has  been 
suggested  that  Darcy’s  Law  is  valid  only  in  those  instances 
where  Re^  <  1  [Ref.  4].  Therefore,  internal  Reynolds  number 
was  calculated  in  all  test  cases  conducted.  In  all  instances 
Re^  was  on  the  order  of  1  or  less. 

Pore  velocity  is  a  critical  factor  in  the  development  of 
this  mathematical  model.  This  air  flow  enhances  internal  con¬ 
vection  heat  transfer,  providing  a  means  for  extinction  and 
eventual  cooling  of  the  graphite  fibers.  The  air  flow  is  also 
a  source  of  oxygen  for  self-sustaining  combustion,  providing 
a  means  for  heat  generation  within  the  porous  structure.  If 
the  heat  transfer,  at  a  point,  due  to  the  air  flow  is  greater 
than  the  heat  generation,  then  cooling  will  occur  at  the  point 
If  heat  generation  is  greater  than  heat  transfer  at  a  point, 
then  combustion  will  proceed  at  this  location. 

Equation  (II.  1]  indicates  that  pore  velocity  (u^)  is  a 
function  of  permeability  (m) .  As  a  result  of  the  critical 
relationship  between  pore  velocity  and  thermal  response,  it 
was  immediately  assumed  that  permeability  would  be  a  critical 
factor  in  the  formulation  of  the  mathematical  model. 

Specific  permeability  is  indicative  of  the  ability  of  a 
porous  medium  to  allow  flow.  It  is  dependent  on  geometry  and 
on  the  properties  of  the  medium  itself.  Several  empirical 
models  have  been  developed  to  obtained  values  for  permeability 
Vatikiotis  selected  a  serial  type  representation  proposed  by 
Scheidegger  [Ref.  5]  in  formulating  this  mathematical  model: 
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(II. 2) 


parameters : 

CD 

Plate  thickness 

CD) 

C2) 

Filament  diameter 

(d) 

CD 

Permeability 

(ra) 

C4) 

Porosity 

(P) 

CD 

Exterior  velocity 

e-u  ■ 

B.  GOVERNING  EQUATIONS 

After  an  expression  for  average  pore  velocity  had  been 
derived,  appropriate  energy  equations  were  developed. 
Vatikiotis  {Ref.  3]  performed  energy  balances  on  both  the 
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porous  solid  and  on  the  air  within  the  graphite  fiber  structure 
in  order  to  develop  the  proper  equations.  The  energy  balances 
utilized  the  control  volume  approach.  This  approach  is  sum¬ 
marized  as  follows: 

Heat  into  +  Heat  Generation  =  Heat  out  +  Increase  in  Internal 
dv  within  dv  of  dv  Energy 

where  dv  =*  dxdydz. 

By  utilizing  the  porous  medium  as  a  control  volume,  the 
fiber  heat  transfer  equation  was  formulated.  Conduction, 
convection,  radiation,  change  in  internal  energy,  and  heat 
generation  rate  per  unit  volume  were  all  considered  in  the 
formulation  of  this  equation. 

In  developing  the  fiber  heat  transfer  equation,  Vatikiotis 
[Ref.  3]  limited  the  model  to  the  kinetic  regime  of  combustion 
and  assumed  that  the  dominant  chemical  reaction  was 
C  +  02  CQ2 

This  idealized  treatment  of  the  actual  reaction  is  a  good 
approximation  during  the  initial  combustion  of  the  graphite 
fibers . 

An  internal  flow  heat  transfer  equation  was  obtained  by 
performing  an  energy  balance  on  a  differential  volume  of  air. 
Heat  conduction  through  the  air,  convection  from  the  fibers 
to  the  air,  and  energy  transport  due  to  air  flow  were  all 
considered  in  the  formulation  of  this  equation. 

In  developing  the  internal  flow  heat  transfer  equation, 
Vatikiotis  [Ref.  3]  utilized  the  perfect  gas  assumption  for 
air  in  order  to  calculate  enthalpy  as  follows: 


dh  =  C  dT 
a  a 


(II. 5) 


where  Ca  is  the  specific  heat  of  air  at  constant  pressure. 

The  effects  of  the  combustion  by-products  on  air  were 
assumed  to  be  negligible  allowing  the  mass  flow  to  be  repre¬ 
sented  as, 


m  =  o  ®u  dydz 
a  a  p 


(II. 6) 


Equations  II. 5  and  II. 6  were  then  utilized  to  determine 
the  energy  transport  due  to  the  air  flow,  m  h. 

In  both  the  fiber  heat  transfer  equation  and  the  internal 
flow  heat  transfer  equation,  thermal  conduction  was  determined 
by  Fourier's  Law  of  heat  conduction  in  the  following  form: 


"cond  ■  -K4  **** 

Similarly,  convection  heat  transfer  was  determined 
equations  by  application  of  the  following  expression: 


(II. 7) 
in  both 


qconv  ^i^g 


Ta)  dA 


(II. S) 


where  T  is  the  temperature  of  the  graphite,  T  is  the  tempera- 
S  a 

ture  of  the  air,  dA  is  the  surface  area  of  the  graphite  in  the 
differential  volume,  and  h^  is  the  internal  convection  heat 
transfer  coefficient.  An  empirical  expression  developed  by 
Yoshida,  Ramaswami,  and  Hougen  [Ref.  6]  was  used  to  determine 


The  final  field  equation  formed  was  the  oxygen  transport 
equation.  In  order  to  formulate  this  equation  a  mass  balance 
was  performed  on  a  differential  volume  of  porous  medium.  The 
mass  balance  was  formulated  in  the  following  manner: 
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Oxygen  into 
dv 


■  Oxygen  out  +  Oxygen  Consumption  +  Oxygen 
of  dv  Accumulation 

Two  mass  transport  mechanisms  were  considered  in  formulat¬ 
ing  this  equation: 

(1)  n\g  -  molecular  diffusion 

(2)  m  -  convective  transport  due  to  mass  flow 

v 

Molecular  diffusion  was  determined  by  Fick’s  Law, 

mB  -  b||  dvdz  (I  I.  9) 

where  B  is  the  diffusivity  of  oxygen  into  the  porous  medium. 

The  convective  transport  term  was  determined  by 

mc  =  oUp<j>dydz  (II. 10) 

Oxygen  consumption  rate  was  determined  by  multiplying 
reaction  rate,  r(T  ,  <j>)  ,  by  the  stoichiometric  ratio  (1/f)  of 
the  reaction.  The  reaction  rate  is  of  the  Arrhenius  type, 


r  -  A(T,  <|0expC-E/RT) 

where  A"1  is  the  characteristic  time  of  the  chemical  reaction, 

E  is  the  activation  energy,  R  is  the  universal  gas  constant, 

T  is  absolute  temperature  and  b  is  the  concentration  of  oxygen. 
The  reaction  rate  was  determined  by  utilizing  an  empirical 
expression  developed  by  Frank-Karaenetski  [Ref.  7]. 

The  three  field  equations  described  above  are  presented 
below: 


Cl)  Fiber  Heat  Transfer  Equation 

3  3T  h . 2  3T 

TxHXj  -Kr)^i  -  J^v)ag  -  Ta)  *  4Hr(Tg,  *)-»sCg5^  UI.ll) 
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Oxvgen  Transoort  Eouation 


d  U  a  T 

Sx<4£  •  ♦>  ’  H 


(11.13] 


Several  of  the  parameters  included  in  equations  II. 11,  11.12, 
and  11.13  are  functions  of  temperature  (i.e.  thermal  conduc¬ 
tivity,  reaction  rate).  This  temperature  dependence  was 
retained  in  the  mathematical  model  and  the  parameters  were 
changed  with  temperature. 


C.  BOUNDARY  CONDITIONS 

In  order  to  make  the  mathematical  model  results  as  mean¬ 
ingful  as  possible,  it  was  desired  that  boundary  conditions 
be  imposed  that  simulated  the  experiments  of  Fontenot  [Ref.  1] . 

Vatikiotis  imposed  the  following  boundary  conditions: 

3T 

CKg  +  Kr)  -jf  -  hx(Tg  -  TJ  +  ae(Tg--  T/)  at  x  *  0  (11.14] 

3T 

(Kg  +  Kr)  -  -  h2  (Tg  -  TJ  -  ere  (T^  -  T^  '*)  at  x  =  L  (II.  IS] 

T  *  T  at  x  *  0  (11.16] 

a  00 

3T  3T 

-at  *  -st  at  x  '  L  f”-1'! 


b||  *  at  x  -0  (11.18] 

||  -  0  at  x  »  L  (11.19] 
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The  above  boundary  conditions  provide  a  reasonable  model  for 
the  one-dimensional  problem.  Boundary  conditions  11.14  and 
11.15  consider  convection  and  radiation  heat  transfer  from  the 
porous  graphite  structure.  Boundary  condition  11.16  sets  the 
air  entering  the  porous  plate  at  ambient  temperature.  Boundary 
conditions  II. IS  and  11.19  are  Dankwert  conditions  [Ref.  8] 
for  flow  through  a  porous  media. 

D.  ADDITIONAL  ASSUMPTIONS 

The  mathematical  model  is  further  based  on  two  additional 

assumptions  not  previously  described. 

(1)  The  temperature  across  each  individual  fiber  is  as¬ 
sumed  to  be  constant.  This  assumption  is  quite  reasonable  due 
to *the  small  diameter  of  the  graphite  fibers. 

(2)  The  temperature  of  the  air  in  each  individual  pore 
is  assumed  to  be  constant. 

These  two  assumptions  should  not  be  interpreted  as  limiting 
the  temperature  from  fiber  to  fiber,  or  from  pore  to  pore. 

The  model  allows  temperatures  to  vary  according  to  the  equations 
described  above. 

E.  METHOD  OF  NUMERICAL  SOLUTION 

The  field  equations  form  a  system  of  three  nonlinear, 
coupled  partial  differential  equations.  A  Galerkin  formulation 
of  the  finite  element  method  is  utilized  to  obtain  solutions 
of  the  field  equations.  Linear  interpolation  fuctions  are  used 
as  the  basis  functions. 
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The  plate  thickness  (L)  is  considered  to  form  a  closed 
domain  (0,  L)  that  is  partitioned  into  (n-1)  elements.  This 
process  results  in  an  n  nodal  point  model.  In  accordance  with 
the  Galerkin  method,  a  system  of  ordinary  differential  equa¬ 
tions  is  finally  obtained. 

The  resulting  system  contains  3n  ordinary  differential 
equations  (graphite  temperature,  air  temperature,  and  oxygen 
concentration  at  each  nodal  point).  These  equations  retain 
the  character  of  the  original  set  of  partial  differential 
equations . 

Initial  tests  were  conducted  using  a  21  nodal  point  model 
(63  o.d.e.).  Several  instances  of  numerical  instability  were 
discovered  for  twenty  one  nodal  point  computer  runs.  Thirty- 
one  nodal  point  and  thirty-five  nodal  point  models  were  tested, 
with  convergence  obtained  within  five  percent.  A  thirty-one 
nodal  point  model  (93  o.d.e.)  was  selected  for  all  subsequent 
test  runs. 
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III.  TEST  PLAN  FORMULATION 


It  was  desired  that  the  thermal  response  of  a  composite 
plate  be  thoroughly  analyzed.  In  order  to  carry  out  this 
objective,  a  test  plan  was  formulated  for  use  with  the  one¬ 
dimensional  combustion  model  described  in  Section  II.  Initial 
efforts  in  formulating  the  test  plan  concerned  selection  of 
input  parameters  for  variation  during  testing. 

A.  EXTERIOR  WIND  VELOCITY  (Uj 

Exterior  air  flow  rate  (U  )  was  immediately  selected  as 
a  test  parameter  because  of  its  obvious  effects  on  average  pore 
velocity  and  because  it  acts  as  a  source  of  oxygen  for  the 
combustion  process.  It  was  felt  that  relative  wind  velocities 
between  ten  and  ninety  knots  would  be  representative  of  the 
velocities  commonly  experienced  on  aircraft  carrier  flight 
decks.  Wind  velocities  of  10,  40,  70,  and  90  knots  were  se¬ 
lected  for  testing. 


B.  PLATE  THICKNESS  CL) 

Plate  thickness  CL)  is  a  critical  parameter  because  it 
also  effects  pore  velocity.  A  consideration  of  Darcy's  Law 
points  out  the  effects  of  plate  thickness  on  pore  velocity. 


u. 


m  Ap_ 
u  L 


(III.l] 


where  £p  is  the  pressure  differential  across  the  plate,  m  is 
the  specific  permeability  of  the  medium,  and  u  is  the  dynamic 
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viscosity.  As  can  be  seen  in  equation  1 1 1 . 1 ,  pore  velocity 
(u0)  decreases  hyperbolically  as  plate  thickness  increases. 

The  effect  of  olate  thickness  on  pore  velocity  was  felt  to  be 
quite  relevant  because  it  could  produce  extremely  high  or  ex¬ 
tremely  low  pore  velocities,  both  of  which  could  result  in 
extinction.  Plate  thicknesses  of  1.0  inch,  0.5  inch,  0.25  inch, 
and  0.125  inch  were  selected  for  testing.  It  was  felt  that 
these  thicknesses  were  typical  of  composite  plates  under  con¬ 
sideration  for  use  as  airframe  material. 

C.  INTERNAL  PLATE  GEOMETRY 

The  pore  velocity  equation  (1 1 . 1 )  indicated  that  both  per¬ 
meability  and  porosity  of  the  porous  medium  were  significant 
factors  in  analyzing  the  thermal  response  of  the  composite 
plate.  Porosity  is  a  measure  of  void  space  and  is  indicative 
of  the  space  available  for  oxygen,  while  permeability  represents 
the  ability  of  a  porous  medium  to  allow  flow.  Equations  II. 2 
and  II. 3  indicated  that  filament  diameter  (d)  and  ply  thickness 
(D)  significantly  effect  both  porosity  and  permeability  (see 
Figure  1).  Filament  diameters  of  .0002  inch,  .0004  inch,  .002 
inch,  and  .004  inch  were  selected  for  testing.  Ply  thicknesses 
of  .0005  inch  and  .005  inch  were  utilized  in  the  testing  program 
These  selections  of  filament  diameter  (d)  and  ply  thickness  (D) 
allowed  for  consideration  not  only  of  the  effects  of  fiber 
diameter,  porosity,  and  permeability,  but  also  the  effects  of 
changing  permeability  while  holding  porosity  constant. 
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It  was  also  desired  to  consider  the  effects  of  changing 
porosity  while  holding  permeability  constant.  In  order  to 
accomplish  this  objective,  specific  filament  diameters  and 
ply  thicknesses  were  selected  such  that  permeability  remained 
the  same  in  each  case.  The  following  selections  accomplished 
this  task: 

1.  d  *  .003  inch;  D  *  .004  inch 

2.  d  3  .006  inch;  D  =  .007  inch 

3.  d  =  .0002  inch;  D  =  .0007  inch 

The  final  test  plan  is  presented  as  Table  II.  This  (4  x  28) 

matrix  of  computer  runs  resulted  in  production  of  112  differ¬ 

ent  cases  for  analysis.  In  addition  to  graphite  temperature 
(T  ),  air  temperature  (T  ),  and  oxygen  concentration,  the 
following  system  parameters  were  calculated  as  functions  of 
time  and  position  for  each  case: 

(1)  Filament  diameter 

(2)  Porosity 

(3)  Permeability 

(4)  Pore  velocity 

(5)  Internal  Reynolds  number 

(6)  Internal  heat  transfer  coefficient 

(7)  Reaction  rate 

(8)  Filament  conductivity 

(9)  Diffusion  coefficient  of  oxygen 

Therefore,  each  of  the  computer  runs  produced  data  concerning 
changes  in  temperature  and  oxygen  concentration  and  changes  in 
system  parameters. 
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TABLE  II:  TEST  PLAN 


PLATE 

FILAMENT 

PLY 

CASE 

CASE 

CASE 

CASE 

THICKNESS 

DIAMETER 

THICKNESS 

A 

B 

C 

D 

SET 

[IN] 

[IN] 

[IN] 

[KNOTS]  [KNOTS]  [KNOTS]  [KNOTS] 

1 

0.5 

.004 

.005 

U  =90 

U  =70 

U  =40 

U  =10 

7 

1.0 

.004 

.005 

00 

U  =90 

00 

U  =70 

U°°=40 

00 

U  =10 

3 

0.25 

.004 

.005 

U  =90 

00 

U  =70 

U  =40 

oo 

U  =10 

4 

0.125 

.004 

.005 

GO 

U  =90 

00 

U  =70 

GO 

U  =40 

CO 

U  =10 

5 

0 . 5 

.002 

.005 

U°°=90 

oo 

1^=70 

ao 

U  =40 

U°°=10 

6 

1.0 

.002 

.005 

U°°=9  0 

U.-70 

U"=40 

u”=10 

*7 

/ 

0.25 

.002 

.005 

U*  =  90 

U  =70 

U  =40 

U.-10 

s 

0.125 

.002 

.005 

U  =90 

Uao=70 

U*=40 

U  =10 

9 

0.5 

.0002 

.0005 

U.-90 

U^70 

U^-40 

U.-IO 

10 

1.0 

.0002 

.ooos 

Uoo=90 

u.-ro 

U  =40 

U  =10 

11 

0.25 

.0002 

.0005 

U  =90 

U  ="0 

U  =40 

uro=10 

12 

0.125 

.0002 

.0005 

u  =70 

U  =40 

u  =io 

13 

0.5 

.0004 

.OOOS 

Uw=90 

U  =70 

U  =40 

U  =10 

14 

1.0 

.0004 

.0005 

U  =90 

uw=70 

U  =40 

U°°=10 

15 

0.25 

.0004 

.0005 

U  =90 

11^=70 

U  =40 

U°°=10 

16 

0.125 

.0004 

.0005 

U>90 

U-?0 

U  =40 

U°°=10 

17 

0.5 

.003 

.004 

1^=90 

00 

U>40 

u“=10 

18 

1.0 

.003 

.004 

U  =90 

U  =70 

U  =40 

U  =10 

19 

0.25 

.003 

.004 

U”=9  0 

U  =70 

U“=40 

00 

U  =10 

20 

0.125 

.003 

.  004 

U  =90 

u“=70 

U“=40 

oo 

U  =10 

21 

0.5 

.006 

.007 

U*=90 

U  =70 

Uoo=40 

u"  =  10 

22 

1.0 

.006 

.007 

U  =90 

U>70 

u  =40 

1^=10 

23 

0.25 

.006 

.007 

U  -90 

U  =70 

U  =40 

Ua5=10 

24 

0.125 

.006 

.007 

U  =90 

U  =70 

CO 

U  =40 

uoo=10 

25 

0.5 

.  0002 

.0007 

U  =90 

U°°=  7  0 

00 

U  =40 

u  =10 

26 

1.0 

.0002 

.0007 

Um=9  0 

u!=^o 

_  00 

U  =40 

oo 

U  =10 

27 

0.25 

.0002 

.0007 

U  =90 

U  =70 

00 

U  =40 

U  *10 

28 

0.125 

.0002 

.0007 

U  =90 

oo 

U*=70 

33 

00 

U  =40 

OO 

oo 

U  =10 

00 

IV.  RESULTS  AMD  OBSERVATIONS 


In  analyzing  each  computer  run,  particular  attention  was 
directed  at  the  combustion-extinction  phenomenon.  Initial 
output  indicated  that  a  particular  combustion  flash  point  did 
not  exist  for  prescribed  initial  conditions.  Instead,  it  was 
surmised  that  a  combustion  continuum  exists  for  a  composite 
plate,  ranging  from  immediate  extinction,  to  immediate  total 
plate  combustion.  This  conjecture  resulted  in  an  additional 
research  objective  of  attempting  to  determine  the  nature  of 
this  continuum.  Placement  of  the  different  computer  runs  in 
the  speculated  combustion  continuum  was  a  necessarily  subjec¬ 
tive  process.  The  following  general  categories  were  utilized: 

(1)  Immediate  extinction 

(.2)  Slight  combustion,  fast  cooling 

(3)  Slight  combustion,  slow  cooling 

(4)  Moderate  combustion 

(.5)  Extensive  combustion 

(6)  Rapid  combustion 

The  effect  of  each  of  the  input  parameters  on  the  combus¬ 
tion-extinction  process  was  studied  extensively.  It  was 
expected  that  by  analyzing  these  effects  for  numerous  computer 
runs,  that  the  combustion  continuum  for  a  particular  plate  and 
prescribed  initial  conditions  would  be  understood,  and  that  the 
effects  of  system  parameters  would  be  discovered.  This 
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information  could  then  be  utilized  to  select  those  design  para¬ 
meters  that  improve  the  flame  resistance  of  composite  materials. 
Improved  flame  resistance  would  increase  the  survivability  of 
a  composite  structure. 

A  graphics  package  was  inserted  in  the  program  to  allow  for 
ease  of  interpretation  of  results.  Graphite  temperature,  air 
temperature,  oxygen  concentration,  reaction  rate,  and  the 
various  system  parameters  listed  above  were  all  included  in 
the  plotting  procedure.  The  plots  are  presented  as  functions 
of  plate  location  for  particular  times. 

The  initial  conditions  prescribed  for  the  112  computer  runs 
are  presented  in  Table  III. 

An  initial  graphite  temperature  of  10S0°F  was  selected  as 
a  result  of  several  preliminary  test  runs.  These  test  runs 
indicated  an  ignition  temperature  in  the  vicinity  of  1050°F 
for  filament  diameters  (d)  and  ply  thicknesses  (D)  of  the  same 
order  of  magnitude  as  those  called  for  by  the  test  plan. 

Initial  air  temperatures  were  selected  to  simulate  a  steady 
increase  in  temperature  from  ambient  (50°F)  at  the  entrance  to 
10S0°F  in  the  plate  interior. 

It  was  observed  that  initial  oxygen  concentration  had 
little  effect  on  the  solution.  The  mathematical  model  imme¬ 
diately  settled  on  the  oxygen  concentration  appropriate  for 
the  particular  test  case,  no  matter  which  initial  oxygen  con¬ 
centration  was  utilized. 


TABLE  III:  INITIAL  CONDITIONS 


n'O'daL 

POINT 

X/L 

Mgi 

~Ain - 

TEMPERATURE 

Oxygen 

CONCENTRATION 

[°F] 

[°F] 

[  Ibm/FT  3 ] 

1 

0.0 

1050.0 

50.0 

.004 

2 

0.03333 

1050.0 

150.0 

.  004 

3 

0.06667 

1050.0 

250.0 

.  004 

4 

0.10000 

1050.0 

350.0 

.004 

5 

0.13333 

1050.0 

450.0 

.004 

6 

0.16667 

1050.0 

5  50.0 

.  004 

7 

0.20000 

1050.0 

6S0.0 

.004 

8 

0 . 23333 

1050.0 

750.0 

.004 

9 

0.26667 

1050.0 

850.0 

.004 

10 

0.30000 

1050.0 

950.0 

.004 

11 

0.33333 

1050.0 

1050.0 

.004 

12 

0.36667 

1050.0 

1050.0 

.004 

13 

0.40000 

1050.0 

1050.0 

.004 

14 

0.43333 

1050.0 

1050.0 

.004 

15 

0.46667 

1050.0 

1050.0 

.004 

16 

0.50000 

1050.0 

1050.0 

.004 

17 

0.53333 

1050.0 

1050.0 

.004 

18 

0.56667 

1050.0 

1050.0 

.004 

19 

0.60000 

1050.0 

1050.0 

.004 

20 

0.63333 

1050.0 

1050.0 

.  004 

21 

0.66667 

1050.0 

1050.0 

.  004 

22 

0.70000 

1050.0 

1050.0 

.004 

23 

0.73333 

1050.0 

1050.0 

.  004 

24 

0.76667 

10S0.0 

1050.0 

.  004 

25 

0.8000Q 

1050.0 

1050.0 

.004 

26 

0.33333 

1050.0 

1050 . 0 

.  004 

27 

0.86667 

1050.0 

1050.0 

.004 

28 

0.90000 

1050.0 

1050.0 

.004 

29 

0.93333 

1050.0 

1050.0 

.004 

3Q 

0.96667 

1050.0 

1050.0 

.004 

31 

1.0000 

1050.0 

1050.0 

.004 

These  initial  conditions  produced  results  that  ranged 
from  immediate  extinction  to  moderate  combustion  for  the  112 
test  runs.  The  influence  of  the  input  parameters  on  the  com¬ 
bustion-extinction  phenomenon  was  quite  pronounced  in  several 
cases,  these  effects  will  now  be  considered. 

A.  EFFECTS  OF  PLATE  THICKNESS 

The  test  plan  provided  four  plate  thicknesses  in  order  to 
explore  the  effects  of  this  parameter  on  thermal  response. 

The  thicknesses  were  as  follows: 


(1) 

0. 

.125 

inch 

(2) 

0. 

.250 

inch 

(3) 

0  . 

.500 

inch 

(4) 

1, 

.000 

inch 

It  was  generally  observed  that  varying  plate  thickness  pro¬ 
duced  two  extinction  regions,  existing  on  both  sides  of  a 
combustion  region.  Slight  to  extensive  combustion  was  observed 
for  the  test  runs  when  in  the  combustion  region.  Th^  two  ex¬ 
tinction  regions  are  the  result  of  two  different  phenomena: 

(1)  'Choking'  of  oxygen  flow  -  pore  velocity  and  oxygen 
concentration  are  depleted  to  such  an  extent  that 
combustion  cannot  be  supported; 

(2)  'Blowout'  -  oxygen  concentration  and  pore  velocity 
are  of  such  magnitude  that  convective  heat  transfer 
overcomes  heat  generation  and  cooling  occurs. 
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As  a  specific  example  of  plate  thickness  effects  consider 
test  sets  1-4.  The  following  parameters  characterized  the  plate 
geometry  for  these  test  sets: 

(1)  Filament  diameter  =  .004  inch 

(2)  Ply  thickness  =  .005  inch 

The  .125  inch  plates  (test  set  4)  experienced  extinction  for  all 
exterior  velocities.  Figure  2  depicts  the  extinction  process 
for  the  0.125  inch  plate,  with  an  exterior  wind  velocity  of  90 
knots.  Graphite  temperatures  are  decreasing  at  all  locations 
and  oxygen  concentration  is  high  at  all  locations,  ranging  from 
.007-. 0j?  Ibm/FT3.  These  thin  plates  are  characterized  by  high 
pore  velocities,  typically  greater  than  10,000  FT/HR.  (see  Figure 
3).  These  high  pore  velocities  and  oxygen  concentrations  indicate 
that  extinction  occurs  as  a  result  of  the  'blowout'  phenomenon. 

The  high  pore  velocities  enhance  convective  heat  transfer  pro¬ 
ducing  the  cooling  effect  depicted  in  Figure  2.  The  high  pore 
velocities  were  expected  for  the  thin  plates  because  pore  velo¬ 
city  is  governed  by  Darcy's  Law.  As  previously  pointed  out,  pore 
velocity  increases  hyperbolically  as  plate  thickness  decreases. 

The  .25  inch  plates  (test  set  3)  experienced  moderate 
combustion  for  all  velocities  tested,  with  the  exception  of 
ten  knots.  Figure  4,  5,  and  6  depict  a  typical  combustion 
sequence  for  plates  of  this  thickness.  As  these  figures 
indicate,  the  combustion  process  initiates  in  a  region  removed 
from  the  plate  entrance,  eventually  moving  to  the  surface  of 
the  plate.  It  appears  that  some  initial  cooling  of  the  graphite 
fibers  is  occurring  in  the  plate  entrance  regions.  The 


situation  is  only  temporary  though,  until  self-sustaining  com¬ 
bustion  initiates  in  the  plate  interior.  The  oxygen  supply 
follows  the  combustion  sequence  with  all  oxygen  being  consumed 
in  the  combustion  region  and  zero  oxygen  concentration  existing 
beyond  this  region  (see  Figure  7).  The  oxygen  concentration 
results  explain  the  movement  of  the  combustion  region  from 
the  plate  interior  to  the  plate  surface.  The  combustion  pro¬ 
cess  is  utilizing  all  the  available  oxygen  in  order  to  sustain 
itself,  but  additional  oxygen  could  be  used  to  enhance  combus¬ 
tion.  Therefore,  the  combustion  region  moves  toward  the 
entrance,  where  the  most  abundant  oxygen  supply  exists. 

Pore  velocities  for  the  0.250  inch  plates  are  greatly 
reduced  from  those  of  the  previously  discussed  thinner  plates. 
Velocities  of  7000  FT/HR  are  typical.  The  increased  thickness 
reduced  the  pore  velocity  to  such  an  extent  that  convecting 
heat  transfer  is  no  longer  removing  heat  at  a  rate  great  enough 
to  provide  cooling,  but  the  velocity  is  of  sufficient  magnitude 
to  provide  enough  oxygen  to  support  moderate  combustion  in  the 
later  half  of  the  plate.  It  is  thought  that  combustion  initiates 
in  the  plate  interior  and  not  on  the  plate  surface  because 
pore  velocity  is  high  enough  to  provide  some  initial  cooling 
at  the  plate  entrance.  It  is  further  thought  that  continued 
reductions  in  pore  velocity  would  result  in  reduction  of  the 
cooled  entrance  region.  A  critical  pore  velocity  could  be 
reached  where  no  entrance  cooling  is  achieved,  but  the  velocity 
is  sufficient  to  support  combustion  throughout  the  entire 


plate.  Pore  velocity  reductions  beyond  this  critical  point 
would  result  in  insufficient  oxygen  supply  to  support  combus¬ 
tion  throughout  the  entire  plate.  Combustion  initiation  in 
the  front  half  of  the  plate  would  characterize  this  situation. 
The  results  obtained  from  the  .5  inch  plates  appear  to  support 
this  speculation. 

The  .5  inch  plates  (test  set  1)  experienced  moderate  com¬ 
bustion  similar  to  the  .25  inch  plates  for  all  velocities 
tested,  with  the  exception  of  ten  knots.  The  combustion  pro¬ 
cess  for  these  plates  differed  from  that  of  the  .25  inch  plates 
in  the  location  of  the  combustion  initiation  region.  As  can 
be  seen  in  Figure  8,  combustion  initiates  in  the  plate  entrance 
region.  Pore  velocities  of  3500  FT/HR  were  typical  for  these 
plates.  It  appears  that  the  reduction  in  pore  velocity  greatly 
reduced  the  cooling  effect  in  the  plate  entrance  regions  noted 
in  the  .25  inch  plates,  but  that  the  velocity  was  still  high 
enough  to  support  combustion  in  this  region. 

The  1.0  inch  plates  (test  set  2)  experienced  slight  local 
combustion  in  the  plate  entrance  region,  with  eventual 
extinction  occurring  for  all  velocities  tested.  In  all  cases, 
pore  velocities  were  not  sufficient  to  support  combustion  for 
an  extended  period  of  time.  As  can  be  seen  in  Figure  9,  oxygen 
concentration  is  rapidly  falling  even  with  an  exterior  wind 
velocity  of  ninety  knots.  This  lack  of  oxygen  within  the 
plate  interior  effectively  chokes  the  combustion  process 
resulting  in  extinction  at  all  velocities  tested. 
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The  effects  of  thickness  on  thermal  response  are  quite 
evident  from  the  above  tests.  Thin  plates  are  desirable  in 
that  extinction  by  'blowout'  is  frequently  achieved,  but  thin 
plates  are  still  susceptible  to  combustion  if  average  pore 
velocity  is  reduced  to  such  an  extent  that  the  thermal  response 
moves  from  the  'blowout'  region  to  the  combustion  region.  Re¬ 
duced  exterior  wind  velocity,  or  a  change  in  internal  plate 
geometry  that  results  in  reduced  permeability  could  both  move 
the  thermal  response  for  a  thin  plate  to  the  combustion  region. 

Thick  plates  are  advantageous  in  that  extremely  low  pore 
velocities  are  frequently  achieved.  This  situation  results  in 
a  'choking'  of  oxygen  flow  that  will  lead  to  extinction.  The 
choking  phenomena  is  not  particularly  desirable  in  that  the 
plate  interior  does  not  cool  as  rapidly  as  when  extinction 
occurs  by  'blowout'.  High  temperatures  remain  in  the  center 
of  the  plate  for  long  periods  of  time.  This  type  of  situation 
probably  occurred  when  'smoldering'  was  observed  in  several  of 
the  tests  conducted  by  Fontenot  [Ref.  2]. 

Thick  plates  are  also  susceptible  to  combustion  if  high 
exterior  velocities  occur  or  if  plate  geometry  is  altered  so 
that  permeability  is  increased. 

In  analyzing  the  tests  for  various  thicknesses,  a  good 
deal  was  learned  about  combustion  initiation.  It  was  learned 
that  the  initiation  region  is  controlled  by  pore  velocity, 
which  is  influenced  by  plate  thickness.  It  was  also  learned 
that  the  size  of  the  combustion  initiation  region  is  quite 


33 


sensitive  to  pore  velocity,  with  a  critical  velocity  existing 
for  a  particular  plate  geometry,  such  that  the  entire  plate 
will  experience  combustion  if  this  velocity  is  present. 

B.  EFFECTS  OF  FILAMENT  DIAMETER 

Altering  fiber  diameter,  while  holding  all  other  system 
parameters  constant,  again  results  in  the  three  characteristic 
thermal  response  regions  previously  described  (choking  -  com¬ 
bustion-blowout)  .  Varying  only  fiber  diameter  results  in  a 
change  in  the  internal  geometry  of  the  plate.  A  decrease  in 
fiber  diameter  (d)  will  effect  an  increase  in  both  porosity 
and  permeability,  resulting  in  more  space  available  for  oxygen 
within  the  porous  graphite  structure  and  less  internal  resis¬ 
tance  to  oxygen  flow.  This  situation  was  explored  in  more 
detail  in  test  sets  5-8.  In  these  test  sets  fiber  diameter 
was  reduced  from  .004  inch  (test  sets  1-4)  to  .002  inch  (test 
sets  5-8) . 

Extinction  by  'blowout'  was  observed  for  all  thicknesses 
tested  at  both  90  knots  and  70  knots  (Cases  5A,  5B,  6A,  6B, 

7A,  7B,  8A,  and  8B) .  Figure  10  (Case  5A)  is  an  example  of  the 
thermal  response  observed  for  these  test  cases.  The  reduction 
in  fiber  diameter  resulted  in  more  oxygen  per  unit  volume  and 
higher  pore  velocities  than  those  observed  for  similar  cases 
in  test  sets  1-4,  This  situation  enhances  convective  heat 
transfer  producing  the  rapid  cooling  depicted  in  Figure  10 
(Case  5A) .  Reduction  of  fiber  diameter  had  a  pronounced  effect 


on  the  0.25  inch  plate  and  the  0.50  inch  plate  at  both  90 
knots  and  70  knots.  When  fiber  diameter  was  .004  inch,  both 
plates  experienced  moderate  combustion.  When  fiber  diameter 
was  reduced  to  .002  inch  the  thermal  response  was  altered  from 
the  combustion  region  to  the  extinction  by  'blowout'  region. 

Moderate  to  extensive  combustion  was  observed  in  the 
following  test  cases: 


(1) 

Case 

5C 

(Uoo=4°> 

L  =  0 . 50 ,  d= . 002 ,  D= . 00  5) 

(2) 

Case 

5D 

O 

rH 

II 

8 

L  =  0 . 50 ,  d= . 00 2 ,  D= . 005) 

(3) 

Case 

6D 

(^=10, 

L=1 . 0 ,  d= . 00 2 ,  D= . 00 5) 

(4) 

Case 

7D 

L=0 . 2S  ,  d= . 00 2 ,  D= . 00 5) 

Figure  11  (Case  5C)  depicts  an  extensive  combustion  situation, 
while  Figure  12  (Case  5D)  presents  a  more  moderate  thermal  re¬ 
sponse.  In  these  situations  it  appears  that  the  higher  porosity 
and  permeability  is  enhancing  the  combustion  process  by  pro¬ 
viding  more  oxygen.  The  external  velocities  and  plate  thick¬ 
nesses  present  in  these  cases  result  in  pore  velocities  that 
support  combustion,  but  are  not  of  sufficient  magnitude  so  as 
to  provide  a  cooling  effect. 

Extinction  was  again  observed  in  the  remaining  test  cases: 

(1)  Case  7C  (U^-40,  L  =  0.25,  d=.002,  D-.005) 

(2)  Case  8C  (Uw-40,  L-0.125,  d».002,  D=.005) 

(3)  Case  8D  (U^-10,  L-0.125,  d-.002,  D-.005) 

In  these  cases  very  high  pore  velocities  (3500  FT/HR  -55,000 
FT/HR)  and  oxygen  concentrations  were  present.  The  extinction 
mechanism  was  once  again  'blowout'.  It  should  be  noted  that 
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the  second  characteristic  extinction  mechanism,  'choking',  due 
to  insufficient  oxygen  to  sustain  combustion,  was  not  observed 
for  these  test  cases. 

It  was  speculated  that  the  'choking'  extinction  region 
could  be  observed  by  conducting  additional  tests  at  greater 
thicknesses,  thereby  decreasing  the  pressure  differential 
across  the  plate.  The  reduced  pressure  differential  would  pro¬ 
duce  lower  average  pore  velocities.  In  order  to  confirm  this 
speculation  additional  tests  were  conducted  at  the  following 
thicknesses : 

(1)  2  inches 

(2)  5  inches 

(3)  8  inches 

Filament  diamter  and  ply  thickness  remained  at  .002  inch  and 
.005  inch  respectively  for  these  tests  and  external  velocity 
was  90  knots. 

The  2  inch  plate  resulted  in  extensive  combustion.  Extinc¬ 
tion  by  'blowout'  had  been  previously  observed  for  all  thick¬ 
nesses  tested  at  90  knots.  This  result  for  the  two  inch  plate 
confirmed  the  expected  combustion  region.  Combustion  initiation 
for  the  two  inch  plate  occurred  in  the  center  region  of  the 
plate,  with  a  cooling  effect  noted  in  the  vicinity  of  the  plate 
entrance.  The  location  of  the  combustion  initiation  region 
indicates  that  the  pore  velocity  is  slightly  higher  than  the 
velocity  necessary  to  support  combustion  throughout  the  entire 
plate . 
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The  five  inch  plate  resulted  in  extensive  combustion  once 
again,  but  combustion  initiation  now  occurred  near  the  plate 
surface.  This  change  in  position  of  the  combustion  initiation 
region  indicates  that  pore  velocity  has  been  reduced  below 
the  critical  point  where  combustion  would  occur  throughout  the 
entire  plate.  These  results  correlate  well  with  the  results 
previously  discussed  in  section  IV  A. 

The  3  inch  plate  resulted  in  an  extinction  situation  as  a 
result  of  a  'choking'  of  oxygen  flow.  As  depicted  in  Figure 
13,  oxygen  concentration  rapidly  decreases.  Slight  combustion 
occurs  near  the  plate  entrance,  but  the  oxygen  flow  is  not 
sufficient  to  support  even  this  slight  combustion  for  an  ex¬ 
tended  period  of  time. 

The  three  additional  thickness  tests  described  above  con¬ 
firm  the  existence  of  two  extinction  regions  and  one  ignition 
region  for  the  more  porous  plate,  as  previously  speculated. 

Tests  sets  5-8  indicate  that  decreasing  filament  diameter 
while  holding  other  parameters  constant  results  in  an  increase 
in  both  permeability  and  porosity.  This  situation  typically 
results  in  either  extensive  combustion  or  extinction  by  'blow¬ 
out',  as  described  above,  for  plates  with  thicknesses  of  less 
than  1  inch.  These  results  would  indicate  that  permeability 
and  porosity  significantly  effect  the  internal  flow 
characteristics  thereby  controlling  the  thermal  response  of 
the  composite  plate.  As  previously  noted,  pore  velocity  is  a 
function  of  both  permeability  (m)  and  porosity  (p) . 


37 


C.  EFFECTS  OF  PERMEABILITY 


In  order  to  analyze  the  effects  of  permeability,  test  sets 
13-16  are  analyzed  in  comparison  with  test  sets  1-4.  The  fol¬ 
lowing  information  contrasts  these  two  test  groups: 

TABLE  IV:  COMPARISON  OF  TESTS  SETS  1-4  WITH  TEST  SETS  13-16 


PARAMETERS 

TEST  SETS  1-4 

TEST  SETS  13-16 

Filament  Diameter 

.004” 

.0004” 

Ply  Thickness 

.005" 

.0005” 

Permeability 

.29  x  10‘9  FT2 

.29  x  10‘ 1 1  FT2 

Porosity 

.497 

.497 

As  can  be  seen  from  the  above  information,  porosity  remains 
constant  while  permeability  decreases  by  a  factor  of  100  for 
test  sets  15-16.  It  was  expected  that  this  significant  reduc¬ 
tion  in  permeability  would  alter  the  thermal  response  of  the 
composite  plate. 

No  combustion  was  observed  for  any  of  the  tests  conducted 
in  test  sets  13-16.  Extremely  low  pore  velocities  were  charac¬ 
teristic  of  all  the  test  cases.  Pore  velocities  ranged  from 
.11  FT/HR  for  the  1.0  inch  plate  with  an  exterior  velocity  of 
10  knots,  to  134  FT/HR  for  a  .125  inch  plate  with  an  exterior 
wind  velocity  of  90  knots.  These  velocities  are  not  of 
sufficient  magnitude  to  support  even  slight  combustion. 
Extinction  occurs  in  all  cases  as  a  result  of  'choking'  of  the 
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oxygen  flow.  The  1  inch,  0.5  inch,  and  0.25  inch  plates,  all 
experienced  zero  oxygen  concentration  within  the  plate  interior. 
The  complete  depletion  of  oxygen  occurred  at  the  locations 
indicated  in  Table  V. 

TABLE  V:  OXYGEN  DEPLETION  LOCATIONS 


PLATE 

THICKNESS 

X/L-ZERO  OXYGEN 

1.0 

inch 

.03333 

0 . 5 

inch 

.03333 

0.25 

inch 

.26667 

Table  V  indicates  that  for  plates  greater  than  0.25  inch 
not  only  is  pore  velocity  insufficient  to  support  combustion, 
but  permeability  and  pressure  differential  across  the  plate 
are  such  that  oxygen  does  not  even  flow  across  the  entire 
plate.  This  is  a  very  desirable  situation  from  the  standpoint 
of  combustion  prevention. 

A  comparison  of  test  sets  1-4  (d=.004  inch,  D=.005  inch) 
with  test  sets  13-16  (d=.0004  inch,  D=.0005  inch)  shows  the 
significance  of  permeability  to  the  combustion  process.  In 
test  sets  1-4,  combustion  was  observed  in  nine  of  the  sixteen 
tests  conducted.  In  test  sets  13-16,  combustion  was  not  ob¬ 
served  in  any  of  the  tests  conducted,  despite  the  fact  that 
porosity  remained  the  same.  These  tests  indicate  that  a  mini¬ 
mum  porosity  is  a  necessary  condition  for  plate  combustion, 
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but  it  is  not  a  sufficient  condition.  Porosity  is  an  indica¬ 
tor  of  space  available  for  oxygen,  while  permeability  is  a 
measure  of  the  ability  of  a  porous  medium  to  allow  flow  (i.e. 
high  permeability  implies  low  resistance  to  flow).  Previous 
tests  have  indicated  the  critical  nature  of  pore  velocity  in 
its  relationship  to  the  combustion  process.  Therefore,  per¬ 
meability  will  have  a  significant  impact  on  the  resulting 
thermal  response  of  the  composite  plate. 

Test  sets  13-16  (d=.0004  inch,  D=.0005  inch)  also  provided 
some  insight  into  the  cooling  process  of  the  composite  plate. 
Figures  14,  15,  16,  and  17  depict  graphite  temperatures  for 
the  1.0  inch,  0.5  inch,  0.25  inch,  and  0.125  inch  plates  re¬ 
spectively.  These  figures  indicate  that  plate  cooling  procedes 
more  quickly  for  the  thinner  plates.  A  consideration  of  the 
pore  velocities  and  the  oxygen  concentrations  that  characterize 
these  test  cases  is  relevant  to  the  cooling  process.  Table  VI 
presents  this  information  for  the  four  plate  thicknesses  tested. 
The  oxygen  concentrations  presented  in  Table  VI  all  occurred 
very  rapidly  after  the  problem  start  time  (t  less  than  1.0 
second) . 

Table  VI  illustrates  pore  velocity  increases  as  plate 
thickness  decreases.  This  situation  occurs  because  as  plate 
thickness  decreases  the  pressure  gradient  across  the  plate 
increases.  The  increased  pressure  gradient  results  in  increased 
pore  velocity.  Despite  the  fact  that  pore  velocity  is  increasing, 
it  is  not  sufficient  to  support  combustion,  even  for  the  0.125 
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inch  plate.  Pore  velocities  of  at  least  500  FT/HR  character¬ 
ized  previous  test  cases  that  experienced  combustion.  As  can 
be  seen  from  Table  VI,  the  pore  velocities  for  the  cases  under 
consideration  are  well  below  this  value. 

TABLE  VI 

PORE  VELOCITY  AND  OXYGEN  CONCENTRATION  (TEST  SETS  15-16) 


PLATE  SIZE 

PORE  VELOCITY 

OXYGEN'  CONCENTRATION 

1.0  inch 

17  FT/HR 

0  at  X/L  =  .03333 

0.5  inch 

35.0  FT/HR 

0  at  X/L  =  .03333 

0.25  inch 

68.0  FT/HR 

0  at  X/L  =  .26667 

0.125  inch 

134.0  FT/HR 

.  005  lbm/FT  3 
(throughout  entire 
plate) 

Even  though  the  pore  velocities  experienced  for  these  test 
cases  were  of  insufficient  magnitude  to  support  combustion,  it 
appears  that  they  were  high  enough  to  provide  some  cooling  to 
the  plate  through  convective  heat  transfer.  In  fact,  the 
cooling  was  quite  significant  for  the  0.25  inch  plate  and  the 
0.125  inch  plate  as  Figures  17  and  18  demonstrate. 

It  appears  that  'choking'  of  oxygen  flow  and  rapid  cooling 
occur  simultaneously  for  thin  plates  with  reduced  permeability 
The  reduced  permeability  provides  the  mechanism  for  preventing 
initial  combustion  by  reducing  pore  velocity  to  such  a  level 
that  combustion  cannot  be  sustained.  The  thin  plate  thickness 
results  in  a  pore  velocity  that  is  sufficiently  high  to  cool 
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the  plate  through  convective  heat  transfer.  These  results  in¬ 
dicate  that  the  composite  plate  designer  can  not  only  design 
for  combustion  prevention,  but  also  for  desired  cooling  rate. 
This  can  be  accomplished  by  not  only  designing  for  a  pore 
velocity  such  that  the  thermal  response  is  in  the  extinction 
by  choking  region,  but  also  designing  the  plate  geometry  such 
that  the  oxygen  flow  is  not  completely  choked  within  the  plate 
interior.  It  is  believed  that  a  threshold  pore  velocity  exists 
for  any  particular  composite  plate,  such  that  combustion  will 
not  be  supported  and  a  maximum  cooling  rate  will  be  achieved. 

In  tests  sets  9-12  (d=.0002,  D=.0005)  the  permeability 

-  1  0  2 

is  slightly  higher  (.625  x  10  FT  )  than  that  for  test  sets 
13-16  (.29  x  10  FT  ) .  This  increase  in  permeability 
resulted  in  higher  pore  velocities  in  all  test  cases.  This 
situation  resulted  in  moving  the  thermal  response  from  the 
extinction  by  'choking'  region  to  the  combustion  region  for 
several  test  cases.  Combustion  was  observed  in  eleven  of  the 
sixteen  tests  conducted  for  these  test  sets. 

In  all  the  test  cases  in  which  combustion  occurred,  ini¬ 
tiation  took  place  in  the  vicinity  of  the  plate  entrance.  It 
was  also  observed  that  the  combustion  process  was  not  of  a 
robust  nature,  with  slight  to  moderate  combustion  observed  in 
each  case.  Figure  18  (case  9C)  presents  the  thermal  response 
of  a  0.5  inch  plate,  with  an  exterior  velocity  of  40  knots, 
that  is  typical  of  the  combustion  process  observed  for  these 
test  sets.  Additionally,  oxygen  concentration  was  rapidly 
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depleted  in  the  vicinity  of  the  plate  entrance  for  all  com¬ 
bustion  situations.  Figure  19  (Case  9C)  depicts  the  oxygen 
depletion  process  observed  in  combustion  situations.  The 
combustion  initiation  region  and  the  rapid  oxygen  depletion 
both  indicate  that  the  combustion  process  is  immediately 
utilizing  all  oxygen  available.  These  results  indicate  that 
the  permeability  is  retarding  the  thermal  response.  A  further 
increase  in  permeability  would  result  in  less  resistance  to 
oxygen  flow,  which  would  produce  pore  velocities  of  greater 
magnitude.  The  higher  pore  velocities  would  move  the  thermal 
response  further  into  the  combustion  region,  resulting  in  a 
more  extensive  combustion  situation. 

It  should  be  noted  that  porosity  was  substantially  higher 
in  test  sets  9-12  (p=.9743)  than  in  test  sets  1-4  (p=.497)  and 
test  sets  13-16  (p=.497).  Previous  tests  indicated  that  the 
thermal  response  was  not  as  sensitive  to  changes  in  porosity 
as  it  was  to  changes  in  permeability.  In  these  tests  though, 
permeability  was  altered  with  porosity.  In  order  to  completely 
investigate  the  effects  of  porosity,  test  sets  17-28  were  con¬ 
ducted.  These  tests  were  designed  to  vary  porosity,  while 
holding  permeability  and  all  other  input  parameters  constant. 

D.  EFFECTS  OF  POROSITY 

Test  sets  17-28  investigated  the  effects  of  varying  poro¬ 
sity,  while  holding  all  other  parameters  constant.  Table  VII 
summarizes  the  parameters  that  characterize  the  various  test 
sets . 
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TABLE  VII:  PARAMETERS  FOR  TEST  SETS  17-28 


TEST  SETS 

FILAMENT 

DIAMETER 

PLY 

THICKNESS 

PERMEABILITY 

17-20 

.003  inch 

.004  inch 

.30  x  10"9FT2 

.  560 

21-24 

.006  inch 

.007  inch 

.30  x  10’9FT2 

.420 

25-28 

.0002  inch 

.0007  inch 

.30  x  10‘9FT2 

.935 

1 .  Test  Sets  21-24  (p=.42Q) 

As  Table  VII  indicates,  porosity  and  permeability  for 
test  sets  21-24  are  the  same  order  of  magnitude  as  those  experi¬ 
enced  in  test  sets  1-4  (m  =  .29  x  10  9FT2,  p  =  .497).  As  ex¬ 
pected,  test  sets  21-24  presented  thermal  responses  very  similar 
to  those  observed  for  test  sets  1-4.  Combustion  was  observed 
in  nine  of  the  sixteen  tests  conducted  for  these  sets,  with 
combustion  initiation  occurring  in  a  region  removed  from  the 
plate  entrance.  It  was  also  observed  that  the  combustion  pro¬ 
cess  was  not  of  a  robust  nature,  with  slight  to  moderate 
combustion  observed  in  each  ignition  case.  Extinction  as  a 
result  of  the  'choking'  of  oxygen  flow  was  observed  for  the 
1.0  inch,  0.5  inch,  and  0.25  inch  plates  when  the  exterior 
velocity  was  ten  knots.  Extinction  as  a  result  of  'blowout' 
was  again  observed  for  the  0.125  inch  plates  at  all  velocities 
tested. 
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2 .  Test  Sets  17-20  (pa.56Q) 

In  test  sets  17-20,  porosity  was  increased  to  .560, 

-  9  2 

while  permeability  remained  at  . 30  x  10  FT  .  This  increase 
in  porosity  (approximately  30  percent)  results  in  more  void 
space  available  for  oxygen.  Combustion  occurred  in  cases 
similar  to  those  that  experienced  combustion  in  test  sets 
21-24.  Extinction  was  again  observed  for  the  1.0  inch,  0.5 
inch,  and  0.25  inch  plates  with  exterior  velocity  equal  to 
ten  knots,  and  for  the  0.125  inch  plates  at  all  velocities. 

The  thermal  response  observed  in  these  test  sets,  when  com¬ 
bustion  took  place,  was  markedly  different  than  that  observed 
for  the  combustion  cases  in  test  sets  21-24.  Extensive  com¬ 
bustion  was  observed  in  several  instances,  especially  for  the 
0.25  inch  and  the  0.50  inch  plates.  Figures  20  (Case  17A) 
and  21  (Case  19A)  are  typical  of  the  thermal  response  observed. 

It  appears  that  the  increase  in  porosity  from  .420  (test 
sets  21-24)  to  .560  (test  sets  17-20)  did  not  induce  combustion 
to  occur  in  any  additional  test  cases,  but  it  did  impact  on 
the  thermal  response.  The  combustion  process  is  of  a  more 
robust  nature.  The  additional  space  available  for  oxygen 
supports  the  combustion  process  once  it  initiates. 

3 .  Test  Sets  25-28  (ps.935) 

In  tests  sets  25-28,  porosity  is  increased  to  .935, 
indicating  a  great  deal  of  void  space  available  for  oxygen. 
Combustion  was  again  observed  in  several  instances.  The  sig- 
nigicant  increase  in  porosity  produced  a  dramatic  change  in 
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the  speed  of  the  observed  thermal  response.  Figure  22  (Case 
25B)  is  indicative  of  the  rapid  thermal  response  observed 
when  combustion  took  place  in  these  high  porosity  tests,  as 
Figure  22  depicts,  graphite  temperatures  in  excess  of  12S0°F 
are  achieved  in  less  than  0.25  seconds.  These  tests  again 
indicate  that  an  increase  in  porosity  does  not  induce  combus¬ 
tion  to  occur  in  any  additional  test  situations,  but  it  does 
impact  on  the  nature  of  the  thermal  response.  High  porosity 
(p  greater  than  0.90)  appears  to  result  in  rapid,  extensive 
combustion  situations  if  permeability  is  such  that  pore  velo¬ 
city  will  support  the  combustion  process  and  allow  combustion 
to  initiate. 

In  summary,  it  is  thought  that  permeability,  and  its 
associated  impact  on  pore  velocity,  determines  whether  the 
thermal  response  will  be  located  in  the  combustion  region  or 
in  one  of  the  two  extinction  regions  that  exist  on  either  side 
of  the  combustion  region  (extinction  by  'choking'  of  oxygen 
flow  or  extinction  by  'blowout').  If  the  thermal  response  is 
located  in  the  combustion  region,  it  is  thought  that  poro¬ 
sity  controls  the  speed  and  robustness  of  the  combustion 
process . 

E.  EFFECTS  OF  EXTERIOR  VELOCITY 

Exterior  velocity  is  the  driving  force  for  the  internal 
oxygen  flow  that  is  required  for  both  the  combustion  and  the 
extinction  processes.  Even  though  exterior  velocity  provides 
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this  driving  force,  neither  extremely  high  nor  extremely  low 
velocities  are  sufficient  to  guarantee  extinction.  Previous 
discussion  has  indicated  that  pore  velocity  is  a  critical  fac¬ 
tor  in  analyzing  the  thermal  response.  Permeability,  plate 
thickness,  and  porosity  are  all  factors  which  can  control  and 
alter  interior  pore  velocity  to  such  an  extent  that  the  thermal 
response  may  be  contrary  to  preconceived  expectations  based 
solely  on  the  knowledge  of  exterior  velocity. 

Exterior  velocity  does  appear  to  effect  the  combustion 
initiation  region.  In  several  test  cases,  with  high  exterior 
velocities  (i.e.  90  knots)  initial  combustion  was  observed  in 
the  latter  half  of  the  plate,  while  in  those  cases  with  low 
exterior  velocities,  initial  combustion  occurred  in  the  proxi¬ 
mity  of  the  plate  surface.  These  results  indicate  that  even 
though  plate  geometry  results  in  a  pore  velocity  that  supports 
combustion,  the  exterior  velocity  is  controlling  where  the 
initial  combustion  will  occur.  High  exterior  velocities  result 
in  oxygen  flow  that  provides  cooling  at  the  plate  entrance. 
Initial  combustion  then  occurs  in  the  latter  plate  regions  after 
the  oxygen  temperature  has  reached  a  sufficient  level  so  that 
heat  generation  overcomes  convective  heat  transfer.  Figure  5 
depicts  combustion  initiation  in  the  latter  plate  regions. 

Low  exterior  velocites  are  not  sufficient  to  produce  the 
initial  cooling  effcts  described  above.  If  the  pore  velocities 
are  sufficient  to  support  combustion,  ignition  will  commence 
in  the  proximity  of  the  plate  surface. 
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In  those  situations  in  which  extinction  occurs  as  a  result 
of  'blowout',  higher  cooling  rates  were  observed  with  higher 
exterior  velocites.  This  result  was  as  expected. 
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V.  CONCLUDING  REMARKS 


The  test  results  indicate  that  permeability,  porosity, 
filament  diameter,  plate  thickness,  and  exterior  wind  velocity 
all  have  significant  effects  on  the  combustion  process.  In 
attempting  to  understand  the  thermal  response  of  a  composite 
plate,  it  must  be  realized  that  two  extinction  regions  exist 
on  either  side  of  a  combustion  region.  Figure  23,  depicts 
this  situation.  For  a  given  temperature,  say  T  =  1050°F,  ex¬ 
tinction  will  result  for  pore  velocities  less  than  u  A  or  for 

P 

pore  velocities  greater  than  u  c.  For  u  A  <  u  <  u  c,  com- 

P  P  P  P 

bustion  will  result  for  all  T  >  1050°F.  Thus  for  a  given 

temperature,  and  fixed  design,  there  is  extinction  at  low  pore 

velocites  (up  <  u  A) ,  followed  by  combustion  for  u  A  <  u  <  u  c, 

P  PPP 

and  again  followed  by  extinction  for  Up>  u^c. 

For  pore  velocities  less  than  UpA,  the  system  is  air  poor, 
and  extinction  occurs  because  heat  generation  is  less  than  the 
heat  transferred  out  of  the  system  as  a  result  of  convection. 

As  Up  is  increased  beyond  UpA,  heat  generation  exceeds  heat 
transferred  and  combustion  with  increasing  system  temperature 
will  occur.  The  minimum  ignition  temperature  occurs  at  up=upg 


Increasing  Up  above  Upg  results  in  increasing  ignition  tempera¬ 
tures.  For  pore  velocities  greater  than  upC>  extinction  will 
occur  once  again  since  heat  transferred  exceeds  heat  generated. 
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At  this  point,  the  system  is  air  rich.  That  is,  there  is  an 
excess  of  oxygen  above  what  is  required  to  sustain  combustion 
and  'blowout'  will  occur. 

It  has  been  established  that  an  ignition  curve  similar  to 
that  discussed  above  for  Figure  23  exists  for  any  composite 
plate.  Average  pore  velocity  is  selected  as  the  abscissa  for 
the  curve  because  of  the  profound  effects  observed  as  a  result 
of  variation  of  this  parameter  during  the  test  runs.  Average 
pore  velocity  is  influenced  by  exterior  wind  velocity,  plate 
thickness,  and  interior  plate  geometry. 

If  a  composite  plate  is  designed  to  retard  combustion  by 
'choking'  of  oxygen  flow,  then  high  exterior  velocities  are 
not  desirable  because  they  may  be  of  such  magnitude  that  the 
thermal  response  is  moved  from  the  'choking'  extinction  region 
to  the  combustion  region.  Similarly,  if  the  designer  intends 
combustion  prevention  to  be  achieved  as  a  result  of  'blowout', 
high  exterior  velocities  would  be  desirable.  In  short,  the 
goals  of  the  composite  material  designer  must  be  understood 
if  combustion  is  to  be  prevented  and  a  desirable  cooling  rate 
of  the  porous  graphite  structure  is  to  be  achieved. 

Additional,  extensive  testing  is  required  in  order  to 
further  quantify  the  thermal  response  of  composite  plates  of 
various  geometries.  Tests  conducted  with  additional  exterior 
wind  velocites  and  under  different  initial  conditions  are  also 
required.  An  exploration  of  the  ignition  process  may  prove 
enlightening.  This  task  could  be  accomplished  by  beginning 


50 


with  an  extinction  situation  and  gradually  increasing 
graphite  temperature  for  each  successive  computer  test 
Combustion  initiation  could  be  observed  along  with  the 
ness  of  the  combustion  process. 
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C.URI:  1.  Idealized  Geometry 
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FIGURE 


CASE  3A 


CASE  3A 
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FIGURE  7.  Oxygen  Concentration  -  Case  3A  (U  =90KTS;  d=.004”-  D=  005 
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CASE  1A 


FIGURE  8.  Graphite  Temperature  -  Case  1A  (Uoo=90KT ;  d=.004";  D=.OOS 
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FIGURE  9.  Oxygen  Concentration  -  Case  2A  (U  =90KTS;  d=.004";  D=.005 
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Graphite  Temperature  -  Case  5A  (Uoo=90KT;  d=.001°;  I)=.005" 
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FIGURE  12.  Graphite  Temperature  -  Case  5D  (Uou=10KT;  d=.002";  D=.OOS" 
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FIGURE  13.  Oxygen  Concentration  (U  =90KT;  d=.Q02";  D=.005" 
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FIGURE  14.  Graphite  Temperature  -  Case  I4A  (U  =90KT;  d=.0004 
D-.0005";  L=1 . 0")  . 
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FIGURE  15.  Graphite  Temperature  -  Case  1 3A  (IM90KT;  d=.00()4' 
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Graphite  Temperature  -  Case 
D= . 0005;  L=0 . 125") . 
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Graphite  Temperature  -  Case  9C  ( IJ  40KT;  <J=.O0O2 
D= .0005" ;  L=0.5"). 
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FIGURE  19.  Oxygen  Concentration  -  Case  9C  (Uco~40KT; 
D=  .000  5” ;  L=0 . 5") . 
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FIGURE  22.  Graphite  Temperature  -  Case  25B  (U  =70KTS;  d=.0002 
D= . 0007"  :  L=0 . 5") . 


TEMPERATURE 


FIGURE  23.  Ignition  Curve  for  a  Composite  Plate 
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